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Abstract: The lowest-order low-energy constants £ and F of chiral pertubation theory 
i can be extracted from lattice data using methods based on the equivalence of random 

matrix theory (RMT) and QCD in the epsilon regime. We discuss how the choice of the 
J> lattice geometry affects such methods. In particular, we show how to minimize systematic 

deviations from RMT by an optimal choice of the lattice geometry in the case of two light 
quark flavors. We illustrate our findings by determining £ and F from lattice configurations 
with two dynamical overlap fermions generated by JLQCD, using two different lattice 
geometries. 
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1 Introduction 

It is well known that QCD in a finite volume V at small quark masses m simplifies as the 
Compton wavelength of the pion, m~ l , becomes large compared to I/ 1 / 4 [1]. In this limit 
the space-time dependence of the low-energy effective theory is suppressed and the theory is 
dominated by the constant mode of the pions. The distribution of the low-lying eigenvalues 
of the Dirac operator can then be calculated in random matrix theory [2], see ref. [3] for a 
review. The low-energy constants (LECs) of chiral perturbation theory are used to map the 
dimensionful quantities of QCD (or the effective theory) to the dimensionless quantities of 
RMT, see, e.g., ref. [4]. Matching lattice data for the low-lying Dirac eigenvalues to RMT 
results then allows for a determination of phenomenologically important LECs. 

The lowest-order LECs are S and F. While £ can be determined rather easily from 
the distribution of the small Dirac eigenvalues, F can be determined only if one includes a 
suitable constant background gauge field [5, 6] such as isospin imaginary chemical potential 
[7, 8] . In the following we discuss the geometry dependence of these methods and show how 
to minimize systematic deviations from RMT by an optimal choice of the lattice geometry. 
We also compare our findings with lattice data of the two-flavor epsilon-regime run of 
JLQCD [9, 10] and extract £ and F from these configurations. 

The paper is structured as follows. In section 2 we briefly review the epsilon expansion 
of chiral perturbation theory at next-to-next-to-leading order (NNLO) which allows for a 
systematic discussion of the geometry dependence of RMT-based methods. In section 3 
we summarize relevant results of RMT for the distribution of the lowest Dirac eigenvalues 
at small imaginary chemical potential. In section 4 we compare the analytic predictions of 
section 2 and 3 to lattice data of JLQCD. We conclude in section 5. 

2 The epsilon expansion at NNLO 

In this section we briefly review the epsilon expansion at NNLO with a small imaginary 
chemical potential ijx, see ref. [11]. In the domain where the Compton wavelength of the 
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pion becomes large compared to V 1 ^, chiral perturbation theory (xPT) can be reordered 
according to the power counting [1] 

V ~ e~ 4 , <9 p ~e, 7r(x)~e, m w ~ e 2 , /i ~ e 2 (2.1) 

with covariant derivative <9 P , pion fields vr(a;), pion mass m^, and chemical potential /u. 1 
The corresponding systematic expansion of xPT is called epsilon expansion. To each order 
in e 2 one can integrate out the space-time dependence and obtain a finite- volume effective 
theory in terms of the constant pion mode. The order in e 2 then translates into the order 
in l/(F 2 \fV). At leading order the finite- volume effective action is given by 

<S L fT° = -\VH Tr(M^U + U l M) - ^VF 2 Tr(CU l CU ) (2.2) 
with constant pion mode 

U = exp[i7r ] , 7T = y J d A x ir(x) , (2.3) 

quark mass matrix M = diag(mi, . . . , mjv>)> and quark chemical potential matrix C = 
diag(/ii, . . . , fJ>Nf)i where is the quark mass and ifj,f is the imaginary chemical potential 
of quark flavor / = 1, . . . , Nf. We find that S^P is identical to the RMT action with nonzero 
chemical potential [8] . Note that the pion decay constant F drops out for vanishing chemical 
potential. At next-to-leading order (NLO) in e 2 the general form of eq. (2.2) remains 
unchanged with E — > Ej^°, F — > F^°, see refs. [12-14] for explicit expressions. In an 
actual lattice simulation we measure effective values E e g and F e g, and we need to include 
finite-volume corrections to recover the infinite-volume values E and F. 

At NNLO and to leading order in the small chemical potential 2 the effective action 
has the form [11] 

S e N ff NLO = -^E e N ff NL ° Tr(Mt[/ + U^M) - ^(F e N ff NLO ) 2 Tr{CU^ l CU ) 

+ T^(VF) 2 Tr(C) [Tr (U Q {M\ C}) + Tr{U^{C, M})] 
+ T 2 E(1/F) 2 Tr({M f , C}U C + {C, MjCU^ 1 

+ {U , Cy^CUoM* + CU {C, U^jMUu 1 ) 
+ T 3 E(yF) 2 TriU^CUoC + C 2 ) Tc{MUq X + M*U ) 
+ T 4 S(yF) 2 Tt{Uq 1 CUoC - C 2 ) Tc{MUq X + M^U ) 
+ T 5 E(VF) 2 Tr([Af f , C]U C + [C, M]CU Q l 

+ [Uo, Cp^CUoM^ + CU [C, U^MUq 1 ) 
+ T 6 (FE) 2 [Tr(MC/ - 1 + M ] U )] 2 + T 7 (yS) 2 [Tr(M?7 ( 7 1 - M ] U Q )] 2 
+ T 8 (FE) 2 [Tr(MCV 1 MCV 1 ) + Tr(M^U M^U )] 

+ H!VF 2 Tt:(C 2 ) + H 2 (VT l ) 2 Tr(M Jf M) +H 3 VF 2 (TrC) 2 (2.4) 

1 In the chiral effective theory a nonzero chemical potential is introduced through the covariant derivative 
of the gauged flavor symmetry (see, e.g., ref. [12]), and therefore the power counting of /x is fixed by the 
power counting of d p . 

2 There are also NNLO terms proportional to V 2 C 4 that have been omitted in (2.4). 
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Figure 1. Finite- volume corrections to £ and F for geometries (a x ) on the left and (b x ) on the 
right with parameters F = 90 MeV, L = 1.71 fm, and m^VV = 1. Taken from ref. [11]. 



with finite- volume effective coupling constants Tj and Hi. The LECs E and F also receive 
further corrections, E — >• E^ LO and -F — )• i^ NL °. The terms in eq. (2.4) that were not 
present in eq. (2.2) cannot be mapped to RMT. These terms are proportional to the Tj 
and Hi . Therefore the magnitude of these coefficients determines the systematic deviations 
from RMT of, e.g., Dirac eigenvalue distributions. The coefficients H\ and do not 
couple to Uq or M and are therefore irrelevant for Dirac eigenvalue distributions (which 
involve derivatives with respect to M in the partially quenched theory). The coefficients 
Ti,n 2 , SS T NLO , and i^ NLO depend on the NLO LECs of xPT and on the geometry of the 
space-time box through finite- volume propagators. Explicit results are given in [11]. 
To be specific we discuss the following lattice geometries from now on, 

(a x ) L = xL, L\ = L 2 = L 3 = L , (2.5a) 
(b x ) L 3 = xL, L = L 1 = L 2 = L, (2.5b) 

where x G {1, 3/2, 2, 3, 4}, and Lj is the extent of the space-time box in direction i (i = 
denotes the temporal direction to which \x couples). In figure 1 we show the finite- volume 
corrections to E and F for the different geometries at NNLO for a set of parameters similar 
to the parameters of the JLQCD two-flavor epsilon-regime run [9, 10]. We note that the 
finite-volume corrections to E are invariant under (a x ) (b x ), while the finite-volume 
corrections to F depend on the choice of geometry. The reason is that the permutation 
symmetry of the four space-time dimensions is broken by the chemical potential, to which 
F couples. For our choice of parameters, geometry (b x ) leads to smaller finite-volume 
corrections to F than geometry (a x ). This was also observed in ref. [12] at NLO. 

We continue our discussion with the finite- volume effective coupling constants Tj and 
%2 that are responsible for the systematic deviations from RMT. The high-energy constant 
H2 should not contribute to low-energy phenomenology (as was shown explicitly for the 
spectral density in [15]), and therefore we do not discuss it further. It is an interesting 
observation [11] that Ti, T2, T3 do not depend on the NLO LECs of xPT and depend on 
the geometry only through a common coefficient 7, i.e., 

Ti,T 2 ,T 3 oc 7 . (2.6) 

The coefficient 7 changes under (a x ) o (b x ), while T4,...,T§ (and %%) are invariant 
under the same exchange [11]. This implies that for nonzero chemical potential a judicious 
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Figure 2. Geometry dependence of systematic deviations from RMT. Taken from ref. [11]. 



choice of geometry is possible which minimizes the systematic deviations from RMT due 
to Ti,T2, T3 at a given volume. Of course, for zero chemical potential (a x ) and (b x ) are 
equivalent, so the choice of geometry has no impact. 

We plot 7 for different geometries in figure 2 for the same set of parameters used in 
figure l. 3 We note that the coefficient 7, and thus a part of the systematic deviations 
from RMT, can be reduced significantly by choosing the geometry (b x ) instead of (a x ) 
for the same value of the asymmetry x. Explicit numbers for the Tj and their impact 
on systematic errors will be discussed in Sec. 4. Note that for large asymmetries x the 
coefficients Tj grow rapidly. Thus for too large values of x the epsilon expansion breaks 
down. This corresponds to the largest individual dimension max(Lj) being significantly 
larger than the Compton wavelength of the pion. 

Note that we can project out a single topological sector by modifying the integration 
domain of the constant pion mode Uq and including a volume-independent determinant 
term in the partition function, see, e.g., refs. [3, 12]. Therefore the discussion of systematic 
deviations from RMT presented in this section is also valid for fixed topological charge. 

3 Random matrix theory 

In this section we summarize some important results of random matrix theory that can be 
used to determine £ and F from fits to Dirac eigenvalue distributions. We consider chiral 
random matrix theory with imaginary chemical potential defined by the partition function 

Z u = I dVdWe~ N ^ wiw+v ^ [] det(£>(4f) + m r f ) , (3.1) 
J f=i 

where m^, . . ., m r N , (ifi^, . . ., i(jT N ) are the masses (imaginary chemical potentials) of the 
sea quarks, see refs. [8, 16], the latter for the case of real chemical potential. The integral is 
over the real and imaginary parts of the elements of the complex N x (N + v) matrices W 
and V with Cartesian integration measure. The random matrix Dirac operator is defined 

3 Since the Ti are of order l/(4ir) 2 we plot 7(47r) 2 , which is then of order 1. 
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W + iffW 




(3.2) 



which has \v\ eigenvalues equal to zero. Therefore v is interpreted as the topological charge. 
Note that rrif and ff are dimensionless quantities. They have to be mapped to physical 
quantities by comparison with the low-energy effective theory of QCD. It was shown in 
ref. [4] that in the limit N — > oo chiral random matrix theory can be mapped to chiral 
perturbation theory using 



rhf = rrifVTi = 2Nm 



H 2 f F 2 V 



2N(f f f 



(3.3) 



where / denotes an arbitrary quark flavor, rrif is the physical quark mass, and fj,f is 
the physical chemical potential. Thus, the low-energy constants £ and F appear in the 
conversion from physical units to dimensionless random matrix units. Note that refs. [4] 
and [8] use a different notation for the dimension of the random matrix Dirac operator. 
The quantities rh j and /t f are often referred to as microscopic scaling quantities due to the 
limit N — > oo. 

The eigenvalue correlation functions for the random matrix model defined by eqs. (3.1) 
and (3.2) in the limit of N — > oo were calculated in ref. [8]. In this section we consider the 
case of Nf = 2 sea quarks with masses rh u and rhd at zero chemical potential. This setup 
corresponds to the two-flavor simulation of JLQCD [9, 10] that is described in more detail 
in section 4. We then compute Dirac eigenvalues x at zero chemical potential and y at 
imaginary chemical potential iS. Note that we could equally well have used a setup with x 
at imaginary chemical potential —i5/2 and y at imaginary chemical potential +i5/2 since 
only the isospin component of the chemical potential is relevant for eigenvalue correlation 
functions [8]. 

We define the two-point correlator 



HZ ~ An(A = 0))S(y - \ m (jl = 5)) 



(3.4) 



where A n = 2N\ r n = \ n VY, and the sum is over all eigenvalues of the random matrix 
Dirac operator at chemical potential fx = and ifi = iS. This correlator allows for a 
discussion of the shift of Dirac eigenvalues due to the imaginary chemical potential id. 
Equation (3.4) is calculated in ref. [8]. The result is given by 



(2) 

/°(l,l)(£>y) = xy det 



J v {im u ) im u J v+ i{im u ) 
J u (irhd) irh d Jv + i{im d ) 



dct 



#11 #12 
#21 #22 



(3.5) 
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where J v is the Bessel function of the first kind, 



*n = det 



12 



21 



22 



det 



det 



det 



l°(x,irh u ) J u (im u ) irh u J v+ i(im u ) 
2°(x,irhd) Ju(im d ) irh d J v+ i{irh d ) 
2°(x,x) Jy{x) xJ u+ i(x) 

l°(x,irh u ) J v {im u ) im u J v+ i(irh u ) 
l°(x,irh d ) J u (irh d ) im d J u+1 (im d ) 
-l-(x,y) e-' 52 / 2 J v {y) e~ p '/ 2 G„(y,6) 

l + (y,im u ) J u (im u ) im u J u+ i(im u 
l + (y,im d ) J u (im d ) im d J v+1 (im d 
l + (y,x) J v (x) xJ u+1 (x) 

l + (y,im u ) J v (im u ) im u J u+1 (im u ) 
l + (y,im d ) J u (im d ) im d J v+1 (im d ) 



(3.6a) 



(3.6b) 



(3.6c) 



(3.6d) 



with 



I°(x,y) 
^(x,y) 
l~(x,y) 



i r 1 

- / dt J v {xVt)J v (yVi) 
1 Jo 



xJ u+1 (x)J u (y) - yJ u+1 (y)J u (x) 



x 2 — y 2 



1 



exp 



dte ±Pt l 2 J v {xSt)J v (yVi), 
xy 



x 2 + y 2 



2 jo 
1 

I 

5 2 1 V 25 2 
Gu(y, 6) = yJ v +i{y) + 5 2 J u {y) , 



S 2 



(3.7a) 

(3.7b) 

(3.7c) 
(3.7d) 



and I v is the modified Bessel function. 

In the limit of small chemical potential 5 2 <C 1 the term proportional to 5~ 2 in Z~ 
dominates. Furthermore, we can perform a large-argument expansion of the Bessel function 
in X~ and ignore all terms of order S 2 , so that 



P?ii){x,y) = H u (x,y,m u ,rh d ) — 



1 



= ex p 



2ttS 2 



(x - yf 
2S 2 



(3.8) 



with 



det 



H u (x,y,m u ,m d ) 



'xy 



l°(y,irh u ) J u (irh u ) im u J u+1 (irh u ) 
2°(y,im d ) J v {im d ) irh d J u+ i(im d ) 
2°(y,x) J v {x) xJ p+ \(x) 



det 



Ju{irhu) im u J u+ i(irh u ) 
Jv{im d ) im d J u+ i(im d ) 



(3.9) 



Note that the prefactor H v is independent of 5. Let us define a probability distribution 
that measures the shift d of the eigenvalues due to the imaginary chemical potential i 5 up 
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to a cutoff x c , 



P(d, x c 



1 



N{x c ) 

H u (d, x c ,m u ,rh d 



(2) 

dx p, x \\(x, x + d) 



exp 




(3.10) 



with 



H v (d, x c ,m u ,rh d ) 
AT(xc) 



1 



N{x c ) 
dd 



dx 



dx H u (x, x + d, m u , rh d ) 

(2) 



p^(x,x + d) 



(3.11) 



The Gaussian factor peaks strongly at d = 0, and thus we can expand H v about d = to 
linear order in d. The constant term in the expansion is fixed by the normalization 



dd P(d, x c 



1 



for <5 2 — > 0. Therefore we have 
P(d,xc) = 



exp 




(l + Cl d + 0(d 2 )), 



(3.12) 
(3.13) 

corresponds to a 



where only c\ depends on x c . Note that to first order in d, P(d, x 
Gaussian distribution with width 5 and center ci<5 2 . 

In section 4 we use a small 5 for the numerical fits. Due to the Gaussian factor, d is 
of order 5, and therefore the contribution of c\ can be neglected (we have confirmed this 
for our numerical results in section 4). We define P(d) to be the leading contribution to 
P(d, x c ) in the limit of small 5, and therefore 

/ d 2 \ 



P(d) 



1 



V2tt5 2 



exp 



25 2 



(3.14) 



This quantity is well-suited to determine 5 and therefore F from a fit to eigenvalue spectra 
obtained in lattice QCD simulations. Note that in the limit of small 5 the distribution does 
not depend on c\. For a related discussion with imaginary isospin chemical potential we 
refer to ref. [17]. 

In refs. [18, 19] the distribution of the lowest Dirac eigenvalue y was calculated analyti- 
cally, and in ref. [20] the calculation was extended to nonzero imaginary chemical potential 
ift. We use the notation of ref. [20]. The distribution of the lowest eigenvalue is given by 



Pi(y) 



a F (0+2) 



(y,o) 



(3.15) 



with gap probability 
2det 



E, 



(0+2) 
0,0 



(y,o) 



Qs(y,rh u ;t = 1) dtQs(y,m u ;t)\ t= i 
Qs(y,m d ;t = 1) d t Qs(y,rh d ;t)\ t =i 



m d I (rh u )h{rh d ) - m u Io(m d )Ii(m u ) 



exp 



(3.16) 
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where 

Qs(y, m;t) = ~ J dr e r ^ 2 I ( V^t m)^^y I x (yft^F)t y) + I (Vt m) . 

(3.17) 

In figure 3 we display P\{y) for different values of fi and m u = rhd- Note that the depen- 
dence on /} and m u is strongly correlated, and therefore it is challenging to use this quantity 
to determine both S and F from a fit to numerical data. Nevertheless, the distribution of 
the lowest eigenvalue is well-suited to determine the scale of y for jl = and therefore E. 



4 Numerical results 



In this section we check the results of section 2 against the epsilon-regime run of JLQCD 
with two dynamical overlap fermions with masses am u = arrid = 0.002 and 16 3 x 32 
lattice points at lattice spacing a = 0.1091(23) fm [9, 10]. For these parameters we 
have m 7r min(Lj) ~ 1, m 7r max(Lj) ~ 2, and rrv^y/V ~ 1.34 (using the GOR value 
m n = 110 MeV, see below). The sea quarks are at zero chemical potential, and topol- 
ogy is fixed to v = 0. We compute the eigenvalues of the valence overlap Dirac operator 
on 460 configurations at zero and nonzero imaginary chemical potential. 4 In this way the 
existing configurations can be used to extract £ and F with low numerical cost. 

We first fit the distribution -Pi (A) of the lowest-lying Dirac eigenvalue at zero chemical 
potential in figure 4 in order to extract the finite- volume effective value a 3 £ e ff = 0.00208(2), 
where we cite the statistical error. This corresponds to the dimensionful value 

S cfr = (231(1)(5) MeV) 3 , (4.1) 

where we cite the statistical error (left) and the error propagated from the uncertainty in 
the lattice spacing (right). The dimensionless value is compatible with a 3 E e fj = 0.00212(6) 
obtained in ref. [9] on the same configurations by a fit to the integrated Dirac eigenvalue 

4 In order to introduce a nonzero imaginary chemical potential ifj, in the overlap Dirac operator we 
multiply the forward (backward) temporal links by a factor of e mM (e _mM ), see [21]. 
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Figure 4. Fit to lowest-lying Dirac eigenvalue distribution -Pi (A) with % 2 /dof = 2.9, a 3 £ e ff = 
0.00208(2). 



distribution. Note that S is renormalization scheme dependent and that we give only the 
values for the lattice renormalization scheme here. Including finite-volume corrections at 
NLO gives an infinite-volume value 

£ = S eff /1.1454 = (221(F)(5) MeV) 3 . (4.2) 

We use only NLO finite- volume corrections here and in the remainder of this paper since 
the NNLO finite-volume corrections and the systematic deviations from RMT are of the 
same order in 1/F 2 y/V. From figure 4 we see that the systematic deviations from RMT 
are significant. As discussed in section 2, they cannot be minimized by a judicious choice 
of lattice geometry. To further reduce the systematic errors in the fit for S, one would 
have to go to a larger volume, which is beyond the scope of this work. Alternatively, one 
could compute -Pi (A) to NNLO, including all non-universal terms. This is a a very difficult 
calculation that nobody has attempted yet. 

Next we fit the shift of the lowest-lying Dirac eigenvalue due to a small imaginary 
chemical potential ifj, in order to extract F as proposed in ref. [7]. As shown in section 
3, RMT predicts a Gaussian distribution with a 2 = fj?F 2 V for the distribution P of the 
difference d between the lowest Dirac eigenvalue at zero and at nonzero imaginary chemical 
potential, see also refs. [7, 8, 17]. In figure 5 we show the resulting fit for geometry (02) 
with finite- volume effective value 

F e ( g 2) = 66(5) (1) MeV, (4.3) 

where we cite the statistical error (left) as well as the error propagated from the uncertainty 
in the lattice spacing (right). We note that the quality of the fit is rather bad (x 2 /dof = 4.2) 
and that this value is not compatible with the result from a fit to meson correlators obtained 
on the same configurations [22], P meS on = 87.3(5.6) MeV. If we include finite- volume 
corrections at NLO we obtain the infinite-volume value 

F (a2) = 50(4)(1) MeV (4.4) 
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Figure 5. Fit to the distribution of Dirac eigenvalue shifts P(d) due to imaginary chemical potential 
a[i = 0.01 with d = dY*V in geometry (a 2 ). The result is given by F cS — 66(5)(1) MeV with 
X 2 /dof = 4.2. We cite the statistical error (left) as well as the error propagated from the uncertainty 
in the lattice spacing (right). 
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Figure 6. Fit to the distribution of Dirac eigenvalue shifts P(d) due to imaginary chemical potential 
a/i = 0.01 with d = dT>V in geometry (b 2 ). The result is given by F cf f = 85(5)(2) MeV with 
X 2 /dof = 0.91. We cite the statistical error (left) as well as the error propagated from the uncertainty 
in the lattice spacing (right). 



so that the agreement is even worse. The bad x 2 /dof = 4.2 suggests that the non-universal 
terms at NNLO, see eq. (2.4) and the subsequent discussion, affect the distribution in a 
non-trivial manner. 

From our discussion in section 2 we learned that we can significantly reduce these 
systematic deviations from RMT by choosing lattice geometry (62) instead of (02). In 
practice this means that we should rotate the lattice by 90 degrees so that we have one 
large spatial dimension instead of a large temporal dimension. In figure 6 we show the 
resulting fit for geometry (62) with good x 2 /dof = 0.91 and 

= 85(5) (2) MeV. (4.5) 
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Including finite-volume corrections at NLO 5 this gives 

= 80(5) (2) MeV, (4.6) 

which agrees within errors with the result from the fit to meson correlators given above. 

The value of x 2 /dof in our fits is stable under variations of the bin size. Specifically, 
if we keep the value of F fixed at the result of eq. (4.5) and vary the bin size from 0.05 to 
0.25, the value of x 2 /dof only changes within about one standard deviation of the x 2 /dof 
distribution, as theoretically expected. Based on the good x 2 /dof in geometry (62), we 
make the hypothesis that in this geometry the non-universal terms in eq. (2.4) are small 
compared to the universal terms. Using the Cramer- von-Mises criterion to estimate the 
goodness of our fit in geometry (62) we obtain a T-value of 0.097, which implies that there 
is no reason to reject our hypothesis. 

Nevertheless, it is instructive to compute the actual values of the Tj which determine 
the size of the systematic errors. One can show that for Nf = 2 the traces of the flavor 
matrices multiplied by T4 and T5 in (2.4) are identical, and therefore (2.4) only depends on 
the combination T4 + T5. Similarly, (2.4) only depends on the combinations T6 + Ts/2 and 
T7 + Tg/2, up to a term proportional to Redet(ML r c {), which for Uq 6 SU(A r j) reduces to 
Re det M and can therefore be neglected for the same reason as for the term proportional 
to H2 - 6 To compute the Tj we use eq. (2.26) of [11] and find 

Ti = -2T 2 = -2T 3 = Pl+ r* = I" ' 033 ge ° metry M ' (4.7a) 

2F A V y 0.0069 for geometry (b 2 ) , 

T 4 + T 5 = -^i = -0.017, (4.7b) 

3pr 1 ir 1 F 

T 6 + T 8 /2 = - 4 \ pA 3 v 4 = -0.0074 , (4.7c) 

T 7 + T 8 /2 = ^|^ = 0.0025, (4.7d) 

see [11] for the notation and the values of PJ 5 . To arrive at the formulas for the Tj in (4.7) 
we used the relation between the three-flavor NLO LECs Li and the two-flavor NLO LECs 
li [23, eq. (11.6)]. To compute the numbers, we took F from (4.6) and lj = I7 = 0.005 [24, 
eq. (19.21)]. For l T 3 and l\ we used the running given in [24, eq. (10.18)] with M = 
and fx = V~ 1 / i and took T3 = 3.13 and I4 = 4.43 [25, table X]. For we used the GOR 
relation, resulting in = y/2mE/F = 110 MeV, where we took £ from (4.2). For the 
combinations of the Tj in (4.7) the dependence on the renormalization scale drops out so 
that the numbers given there are scale independent. 

Note that in geometry (62) all coefficients Tj are quite small, which is consistent with 
our hypothesis that in this geometry the systematic errors due to the non-universal terms 
are under control. In geometry (02) the Ti ; 2,3 are larger by a factor of ~ 5. These terms 



5 If we use finite-volume corrections at NNLO the value is further reduced by 2%. 

6 For fixed topology we have Uq £ U(7V/), and the contribution proportional to Tg Redet(Mt/g) in (2.4) 
cannot be neglected. However, for two flavors Ts does not depend on the geometry [11] and is proportional 
to Lg, which in turn is a combination of two- flavor LECs and HECs that are all of order 1/(4tt) 2 [24]. 
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appear to be the reason for the systematic deviations observed in figure 5. As in the case of 
S, the systematic errors can be reduced either by going to a larger volume, or by computing 
the spectral correlation functions corresponding to the effective action of eq. (2.4) including 
all non-universal terms. Again, this is beyond the scope of this work. 

To estimate the systematic errors on our determination of X and F, we recall that the 
NNLO finite- volume corrections and the systematic deviations from RMT are of the same 
order in 1/F 2 \/V, see the comment after (4.2). The systematic errors due to the NNLO 
finite- volume corrections can be estimated by comparing the NLO and NNLO values given 
in [11, table 2], and we obtain ~ 2% for F and ~ 5% for E. Assuming that for geometry 
(62) the systematic errors due to non-universal contributions in eq. (2.4) are roughly of the 
same numerical size, we assign a total systematic error of ~ 4% to F and of ~ 10% to S. 

5 Conclusions 

We have shown that the geometry dependence of the Dirac eigenvalue distributions at 
nonzero chemical potential strongly influences the determination of F from RMT fits. 
Making a judicious choice of the lattice geometry (in our case, by exchanging the temporal 
axis with one of the spatial axes), this dependence can be significantly reduced such that 
the systematic error on F is kept under control. This makes the RMT-based method 
proposed in ref. [7] a useful alternative to other lattice methods. 

Our final results for £ and F obtained from the two-flavor epsilon-regime run of JLQCD 
are given by 

^ m (2 GeV) = ( 2 GeV ) S = (230(1) (14) MeV) 3 , F = 80(5) (5) MeV , (5.1) 

where both values include finite- volume corrections at NLO and Z^(2 GeV) = 1.14(2) 
[9, 10]. The left bracket gives the statistical error, the right bracket gives the combined 
systematic error (including the systematic error of the conversion to MS and the uncertainty 
in the lattice spacing). The individual systematic errors were added linearly. 

To reduce the systematic errors due to the finite volume, it would be beneficial to repeat 
this study at a larger simulation volume. Also, to eliminate the contamination caused by the 
non-universal terms, one could attempt, within the framework of the finite-volume effective 
theory of ref. [11], to calculate Dirac eigenvalue distributions beyond RMT including the 
systematic deviations at NNLO in the epsilon expansion. Work in this direction is in 
progress. 
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